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^ Abstract: We use one-loop SU(2) L x SU(2) R chiral perturbation theory (SU(2) ChPT) 

to study the behaviour of the form-factors for semileptonic K — > tt decays with the pion 
mass at q 2 = and at <?n iax = (rriK — m-iv) 2 , where q is the momentum transfer. At q 2 = 0, 

' the final-state pion has an energy of approximately m^/2 (for txik S> m n ) and so is not 

soft, nevertheless it is possible to compute the chiral logarithms, i.e. the corrections of 
j— i ■ 0(m 2 log(m 2 )). We envisage that our results at q 2 = will be useful in extrapolating 

. lattice QCD results to physical masses. A consequence of the Callan-Treiman relation is 

that in the SU(2) chiral limit (m u = = 0), the scalar form factor /° at g^ ax is equal to 

\ f( K i I f, the ratio of the kaon and pion leptonic decay constants in the chiral limit. Lattice 

results for the scalar form factor at Q^ax are obtained with excellent precision, but at the 
masses at which the simulations are performed the results are about 25% below and 

' are increasing only very slowly. We investigate the chiral behaviour of f°(q 2 aSL ^) an d hrid 

large corrections which provide a semi-quantitative explanation of the difference between 
the lattice results and f^/f. We stress the generality of the relation /p^^f/max) = / If 
in the SU(2) chiral limit, where P = K,D or B and briefly comment on the potential value 
of using this theorem in obtaining physical results from lattice simulations. 
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1. Introduction 

One of the most precise methods to extract the V us element of the CKM-matrix is to 
use K — ► nivp semileptonic decays (-fQ 3 decays), where i is an electron or a muon. The 
combination |T4 S / + (0)| can be determined from the experimental rate 

r K ^iu e = C 2 K ^^-IS EW (1 + 2A SU(2) +2A EM ) |K S | 2 |/ + (0)| 2 , (1.1) 

where / is a phase space integral which can be evaluated from the experimentally deter- 
mined shape of the form-factors and A SU ( 2 ), Aem and Sew contain the calculable cor- 
rections due to isospin breaking, electromagnetic and short-distance electroweak effects 
respectively. C\ = 1/2 (1) is the Clebsch-Gordan coefficient for the neutral (charged) kaon 
decay and / + (0) is the form factor defined from 

(tt(p«)\u^s\K( Pk )) = ( P K+P,)^f + (q 2 ) + (PK-Pn)^r(q 2 ) (1-2) 



{PK+P.r-^ 7 ^^ 



where q is the momentum transfer q = pk — Pw At q 2 = we have / + (0) = /°(0) . The 
Particle Data Group(2008) [1] quotes 

|K S / + (0)| = 0.21668(45), (1.3) 
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so that in order to obtain \V US \ we need to determine / + (0). In the last four years, fol- 
lowing ref. [2], lattice QCD calculations of / + (0) have been undertaken using dynamical 
simulations with Nf = 2 or Nf = 2 + 1 flavours of sea quarks [3-7], thus enabling the 
extraction of V us . 

In this paper we investigate the behaviour of decay amplitudes with the masses of 
the u and d quarks. This is an interesting problem in itself, but the immediate motivation is 
the need to extrapolate the lattice results for / + (0), obtained with u and d quarks heavier 
than the physical ones (typically with pions with masses in the range m n > 300 MeV), to 
their physical values. 

The mass dependence of the form factors is studied below using SU(2)l x SU(2)# chiral 
perturbation theory (ChPT) at next-to-leading order (NLO) 1 . Conventionally, following 
the seminal paper of Gasser and Leutwyler [8], it has been SU(3) ChPT which has been 
applied to the study of decay amplitudes. However, following the study of the quark 
mass dependence of physical quantities computed in a lattice simulation using Domain 
Wall Fermions [9], together with our colleagues from the RBC and UKQCD collaborations 
we concluded that it may be better to use SU(2) ChPT, at least for some quantities. 
This conclusion is primarily based on the large one-loop effects in SU(3) ChPT found for 
the leptonic decay constant of 'pions' with masses in the range in which the simulations 
were performed. Note also that the strange quark mass (m s ) in lattice simulations can be 
chosen to be at its physical value and so ChPT is not needed to perform the corresponding 
extrapolation (although, since the bare strange quark mass is chosen before the simulation 
is undertaken, in practice there may have to be a small extrapolation to correct for the 
difference between the m s used in the simulation and its physical value; SU(3) ChPT may 
provide useful guidance for this). Of course, in using SU(2) rather than SU(3) ChPT we 
sacrifice some symmetry and therefore some information. 

In ref. [9], together with the RBC and UKQCD collaborations, we developed and used 
SU(2) ChPT for kaon physics; in particular we studied the dependence on the pion mass 
of the mass of the kaon mx, the leptonic decay constant fx and the 5^-parameter which 
contains the non-perturbative QCD effects in K°-K° mixing. We recall the main features 
of the formalism in section 2. An important difference between SU(3) and SU(2) ChPT 
is that with SU(2) ChPT, powers of m^/A 2 , where fhx is the mass of the kaon in the 
limit m u = rrid = and A x is the scale of chiral symmetry breaking, are absorbed into the 
low-energy constants (LECs). In SU(3) ChPT at n-loop order on the other hand, there 
remain errors of 0((rhx / 'A x ) 2(n+1 )) I. The corresponding uncertainties in SU(2) ChPT are 
of 0((m w /A x ) 2 ( n+1 )) and 0((m Jm^) 2 ( n+1 )). 

In this paper we study two aspects of the chiral behaviour of the form factors: 

1. The behaviour of / + (0) = /°(0) with m u and m^. In order to determine V us 
from eq. (1.3) we need the form factor at q 2 = for physical values of the quark 
masses. The result presented in eq. (3.9) below represents the behaviour of the form 

1 For compactness of notation in the remainder of this paper we refer to SU(w)l x SU(n)ij ChPT (for 
n = 2 or 3) as SU(n) ChPT. 
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[MeV] 


?Lx [GeV 2 ] 


f°(q 2 ) 


671(11) 


0.00235(4) 


1.00029(6) 


556(9) 


0.01252(20) 


1.00192(34) 


416(7) 


0.03524(62) 


1.00887(89) 


329(5) 


0.06070(107) 


1.02143(132) 



Table 1: Results for f (qf nax ) from ref. [3] at four values of the quark masses, corresponding to 
the pion masses given in the first column. 



factor with the pion mass at NLO (one-loop order) and can be used to extrapolate 
the lattice results obtained at larger values of m u = to the physical point. 

In order to derive eq. (3.9) one has to overcome a subtlety. At q 2 = 0, 

2^ • PK = m 2 K + ml = fh\ + 0{m\ ) , 

so that E T , the energy of the pion in the rest frame of the kaon, is approximately 
equal to m^/2 and is not small (i.e. it is not of Ofa^)) for m? K 3> m 2 . Since 
SU(2) ChPT is an expansion in powers of masses and momenta of the pions, the 
fact that the external pion in K — > it semileptonic decays is hard complicates this 
power counting. Nevertheless, by integrating by parts, we show in section 3 that an 
expansion in small masses and momenta of 0(m n ) is possible and results in eq. (3.9). 
This is possible because the chiral logarithms arise from soft regions of phase-space 
for the internal pions. 

2. The behaviour of /°(?max) with m u and m^. The maximum physical value 
of q 2 is (mjf — m n ) 2 , corresponding to the pion and kaon both at rest. Using the 
double ratio techniques proposed in ref. [2], f°(q 2 nSLX ) is evaluated with remarkable 
precision in lattice simulations. For illustration we reproduce in table 1 the results 
from the RBC and UKQCD collaborations' simulation on a 24 3 spatial lattice [3]. 
The point which we particularly wish to underline here is that in the SU(2) chiral 
limit (m u = = 0), the Callan-Treiman relation [10] implies that 

AO (1-4) 

where and / are the kaon and pion leptonic decay constants in the SU(2) chiral 
limit. The Callan-Treiman relation was derived for the unphysical value of q 2 = 
m 2 K — m 2 , nevertheless, in the chiral limit it also holds for q^ax- We shall show 
however, that the corrections to the relation are of 0{m- K ) and not the standard 
ChPT corrections of 0(m 2 ). The ratio of physical decay constants is fx/ fir — 1-2 
and in the chiral limit it is a little larger, e.g. the lattice study of ref. [9] finds a 
ratio f^/f — 1.26 2 . In ref. [3] the entries in table 1 were obtained with a strange 

2 Enno Scholz private communication. This particular result is not quoted directly in [9]. 
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quark mass which is a little larger than the physical one and the corresponding value 
of f( K ^/f is about 1.28. This is the value which we use in the numerical estimates 
below, together with / ~ 115 MeV, which is the central value found in [9]. We restrict 
the comparison of ChPT with table 1 to the entries with = 329 and 429 MeV, 
since our experience from ref. [9] is that one-loop ChPT is less reliable at the heavier 
masses. The values of /°(<?max) i n table 1 are equal to 1 within 2% or so, and although 
they are increasing as the quark masses are reduced, the observed increase is very slow 
indeed. As m w decreases from 670 to 330 MeV, /°(^ ax ) increases only from 1.00 to 
1.02 which is still a long way from the expected value of about 1.28 in the chiral limit. 
We investigate the chiral behaviour of / (o , I 2 nax ) up to one-loop order in section 4 and 
find that the chiral logarithm has a large coefficient but the wrong sign to account for 
the extrapolation to f^/f- The coefficient of the linear term in m n , which is not 
calculable in SU(2) ChPT, can be estimated by converting the SU(3) results of ref. [8] 
to the SU(2) theory. We find that it is large with the correct sign but predicts too large 
a ratio between /°(^ ax ) in the SU(2) chiral limit and at the masses where lattice 
simulations are performed. We also study the full prediction from SU(3) ChPT, 
which reproduces qualitatively (and semi-quantitatively) the observed behaviour. 

The plan for the remainder of this paper is as follows. In the following section we briefly 
recall some of the main features of SU(2) ChPT for kaon physics. Within this context, we 
also derive eq. (1.4). Sections 3 and 4 contain the studies of the chiral behaviour of the form 
factors at q 2 = and g^ax respectively. Our calculations have some overlap with those of 
semileptonic decays of S-mesons [11-17] and we discuss the similarities and differences in 
section 5. 

2. SU(2) Chiral Perturbation Theory for Kaons 

We start by briefly summarising the formalism introduced in section II. B of ref. [9] which 
we apply in the following sections to decays. We write the pion matrix, the quark 
mass matrix and the kaon fields in the form: 

We work in the isospin limit so that mi represents m u = m&. The pion matrices £ and £ 
are defined in the standard way: 

£ = exp(^//) and £ = £ 2 , (2.2) 

where / is the pion decay constant in the SU(2) chiral limit, m u = m^ = 0. As with all 
LECs in SU(2) ChPT, / depends on m s , the mass of the strange quark. Throughout this 
paper we define the pion and kaon decay constants using a normalization in which the 
physical value for the pion is f n ~ 131 MeV. 

We need to construct the chiral Lagrangian and operators which transform in a speci- 
fied way under SU(2) chiral transformations out of the fields in eqs. (2.1) and (2.2). Under 
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global left and right handed transformations, L and R respectively, these fields transform 
as follows: 

£ -► L£f7 f = L^R 1 ", S -► LSi? f and K ^UK, (2.3) 

where C/ is a function of L, i? and the meson fields which reduces to a global vector 
transformation when L = R. From the transformations in eq. (2.3) we construct operators 
with the required flavour and chiral quantum numbers. 
The pion Lagrangian at lowest order is well known: 

= ^ tr {fyE^Et} + Z!^ tr { M t£ + MSt} , (2.4) 
8 4 

where -B is the standard lowest order LEC and to this order m 2 = 2Bmi. For the interac- 
tions of kaons, which are not considered soft in the SU(2) ChPT formalism, with soft pions 
the chiral Lagrangian has been introduced by Roessl [18] and at lowest order is given by 

L^ K = D^tf D»K - m 2 K K^K , (2.5) 

where the covariant derivative is constructed using the vector field V^, 

V* = \ + ^) - UVprf + Udntf, (2.6) 

and is defined by 

D^K = d^K + V^K UD^K . (2.7) 
In the following it will be necessary also to introduce the pion axial vector field defined 

by 

A^= l - (t}d^ - - UAprf . (2.8) 

When constructing Feynman diagrams from the irK Lagrangian and the effective theory 
local operators, we expand the vector and axial fields in terms of the pion fields, 

y^ = Jp[^d^} + --- and A lt = -jd lt <f> + --- (2.9) 

so that the first term in the expansion of the vector field contains two pions and that for 
the axial field starts with a single pion. 

Similar calculations for B — > n and B — > K semileptonic decays using heavy meson 
ChPT were undertaken in ref. [15, 16] and we will discuss the similarities and differences 
with kaon decays in more detail in section 5. Here we simply point out that in the heavy 
meson ChPT, the limit m& — > oo is taken before performing the chiral expansion. The 
resulting spin symmetry implies that the B* vector meson is degenerate with the pseu- 
doscalar B, and so the B*Btt interactions (where the pion is soft), and hence diagrams 
containing B* propagators, must also be included. In our kaon case, the K* — K mass 
splitting is considered to be of O(Aqcd) and so the corresponding diagrams are absent. 
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2.1 Lowest Order AS = 1 vector and axial currents 

We end this section by discussing the lowest order AS = 1 vector and axial currents in the 
effective theory. As already mentioned in the Introduction, it will not be enough to consider 
only the lowest order contributions and we will have to extend the present discussion in 
the following two sections. 

The left handed QCD AS = 1 current is 

Jjl = QlI^sl = q^ {1 ~ 2 l5) s , (2.10) 

where q = u or d. It is convenient to promote q to be a 2-component vector with components 
u and d and to introduce a 2-component constant spurion vector h in order to be able to 
project u and d as required; specifically we write the left-handed current as 

qh^Q-^s. (2.11) 

The current in eq. (2.11) would be invariant under SU(2)l transformations if h transformed 
as h — ► Lh. We now construct the form of the left-handed current in the effective theory. 
This is a linear combination of all operators which are linear in h and which would be 
invariant under SU(2)l transformations if h transformed as above. At lowest order in the 
chiral expansion we identify two possible independent terms and, following the notation of 
ref. [9], we write the left-handed current as 

= -L A1 (D^^h + iL A2 K^A^h , (2.12) 

where Lai and La2 are LECs and is the pion axial current defined in eq. (2.8). Note 
that since = —iD^ (and also A^ = iD^), no new independent operator is obtained 
by replacing A^ by the covariant derivative in the second term on the right-hand side of 
eq. (2.12). 

For the right-handed current, we take the transformation on h to be h — > Rh and 
obtain two possible operators at lowest order, 

(D^K)^h and K^A^h, (2.13) 

which transform as ^qhj^l + 75) s. Noting that parity transformations, under which 

K^-K, £-£L, A^-A^ (2.14) 

transform into (where is the right-handed current), so that the same LECs, Lai 
and La2 appear in the right-handed current, 

J* = L A i(D^K)^h + iL A2 K^A^h . (2.15) 

In some applications, the La2 term can be considered to be sub-leading since the derivative 
is on the pion rather than the kaon field. The vector (J M ) and axial-vector (jjj) currents 
can now readily be determined: 

J M = J* + Jf: = LAiiD^tt - £})h + iL A2 K^A^ + g)h (2.16) 
4 = J* ~ Jjl = L A i{D„K)\i + $)h + iL A2 K^A.H - $)h . (2.17) 
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The LEC Lai appears in both the vector and axial-vector currents and we will see in 
section 4 that it is this feature which allows us (in the SU(2) chiral limit) to relate the 
K — ► vacuum matrix element of the axial-vector current and the K — > it matrix element 
of the vector current and hence to derive eq. (1.4). Evaluating the K — > vacuum matrix 
element in the chiral limit immediately shows us that 

= 2L A1 . (2.18) 

The symmetry arguments used here apply also to other flavours, so that eq.(1.4) can 
be generalised to D and B mesons and this is briefly discussed in section 5 below. 

3. K i3 form factors at q 2 = 

As discussed in the Introduction, in order to study the chiral behaviour of the form factor 
f° at q 2 = we have to deal with the fact that in this case 2px • Pw — rn 2 K and hence 
we cannot neglect operators with an arbitrary numbers of derivatives on the external pion 
field. This situation is reminiscent of the light-cone dominated process of deep-inelastic 
scattering. To illustrate the point consider the matrix element 

( Tr( Pn ) | (D v D li K?D''(S - g)h | K( PK ) > (3.1) 

at q 2 = 0. In spite of the additional derivative acting on the pion field relative to the 
first term in in eq. (2.16), the matrix element in (3.1) does give a leading contribution 
in the chiral expansion since by inspection we see that there is a contribution of P k ■ 
p w times the matrix element of (D^Ky^ — £')h. Nevertheless, as we now show, the 
leading contribution is simply proportional to the matrix element of {D^K)^{^ — £')h 
(with a constant of proportionality which depends on m s but not on m u ^) and so the 
chiral logarithms are the same and the number of LECs remains the same. To see this, 
note that at q 2 = 0, 

= d 2 {ir( Pn )\(D^K)\ti-^)h\K(p K )) 
= { ir( Pn ) | (D 2 D^K)\i - $)h + 2{D v D^D v {i - $)h 

+ (D fl K)^D 2 (Z-^)h\K(p K )) (3.2) 

so that 

( ir( Pn ) | {D u D^D%i - tf)h | K( PK ) ) = (3.3) 
-\ {(4) I (D 2 D„K)\t - ?)h | K{p K ) ) + ( n(p w ) | (D.K^D 2 ^ - ?)h \ K{p K ) )} . 

Before discussing the chiral behaviour of the operators on the right-hand side of eq. (3.3) we 
clarify our power counting. The external pion is hard in the sense that 2p n • px — m 2 K and 
so we need to keep p n - P k to any power. This is the reason the matrix element in eq. (3.1) 
is of leading order. We accept that the corrections of 0(m 2 ) are multiplied by an unknown 
constant and so we do not attempt to calculate such terms. We do however, calculate the 
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chiral logarithms, i.e. the corrections of 0(m 2 log (m 2 )) and in order to evaluate these we 
can treat the internal pion momenta as being soft, i.e. of 0(m n ). 

The operator in the second term on the right hand side of eq. (3.3) contains the 
insertion D 2 (£ — £t) . This leads to a contribution which is suppressed by a factor of m 2 , 
with no chiral logarithm proportional to m^log(m^). Thus, up to the order to which we 
are working, we only need to consider the first term on the right-hand side of eq. (3.3), 
where we can replace (D 2 Dp L K)^ by (D^D 2 K)^ = —fh 2 K (D^Ky up to terms which are 
suppressed by m 2 . Note that the commutator [D^,D U ] contains two derivatives acting 
on two different pion fields, at least one of which must be Wick contracted to give a soft 
internal propagator. This leads to a suppression of 0(m 2 ) without chiral logarithms and 
we arrive at the useful result that, up to corrections of Ofa 2 ,) (without chiral logarithms): 

( 7T( Pn ) | {D„D lt K?D''(Z - ?)h | K(p K ) > = ( 7r( P7T ) I {D,K)\H - g)h \ K{p K ) > . 

(3.4) 

Thus in order to include the contribution of the matrix element (3.1) to the K — > ir form 
factor /°(0), including the one-loop chiral logarithms, it is sufficient to replace the LEC 
Lai in the definition of the vector current in eq. (2.16) by an unknown coefficient which 
depends on m s but not on the light-quark masses. 

The discussion of the matrix element in (3.1) presented explicitly above can be gen- 
eralised to other operators. Leading-order operators can have any number of covariant 
derivatives on the external pion field. If the Lorentz index of a covariant derivative acting 
on the external pion field is contracted with another derivative on the external pion then 
we obtain a non-leading correction of 0(m 2 ). Similarly if it is contracted with a derivative 
on a pion in an internal loop we also obtain a similar suppression. Finally if it is con- 
tracted with a derivative on the kaon field then we can reduce it to an operator which is 
proportional to the leading operator by integrating by parts as above. Note also that the 
kaon mass-squared, m 2 K , has no chiral logarithms of the form m 2 log(m 2 ) and so no chiral 
logarithms are introduced by using the equations of motion. 

From this discussion we see that to leading order at q 2 = we have 

(ttCPtt) \q^(l -15)s\K(j>k)) = (3.5) 
( n(p n ) | L A i{D^K)\i - $)h + L A2 K^A^ + g)h \ K(p K ) ) 

where we recall that Lai and L A 2 are unknown constants which cannot be obtained from 
Lai and L A 2 alone. They depend on m s but not on the light-quark masses and hence we 
treat them as LECs, noting however that they are only relevant for the case q 2 = 0. As 
a result of the fact that the matrix element at q 2 = is written in terms of Lai and La2 
rather than Lai and L A 2, we lose the connection to f^/f in this case. 

3.1 The Chiral Logarithms 

The tree level contribution to the matrix element in eq. (3.5) is 

<XPtt) I <77 M (1 - j 5 )s\K(p K ) ) = —j^-pKv + —f^PTT^- (3.6) 
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Pk Ptv Pk Pit Pk Piv Pk Piv 

(a) (b) (c) (d) 

Figure 1: Diagrams contributing to the K — > tt matrix elements at tree level (diagram (a)) and at 
one- loop level (diagrams (b), (c) and (d)). The grey circle represents the insertion of the K — > tt 
vector current and the grey box the insertion of the KKtttt vertex (diagram (b)) or the four-pion 
vertex (diagram (d)) from the strong Lagrangian. 



This contribution is represented diagrammatically in fig. 1(a). 

In order to obtain the chiral logarithms at one-loop order we need to evaluate the 
diagrams in fig. 1(b), (c) and (d), where fig. 1(d) represents the contribution to the pion's 
wave-function renormalization. There is no one-loop chiral logarithm contributing to the 
kaon's wave function renormalization which can be deduced from the structure of the 
D^KD^K term in the strong Lagrangian L^h in (2.5). The KKtttt vertex arises when 
one keeps a partial derivative from one of the factors and the current from the 
other. From eq. (2.9) we see that the expansion of starts with two pion fields, on one 
of which there is a single derivative. This derivative corresponds to a single momentum in 
the numerator of the tadpole loop and hence the momentum integration is odd and gives 
zero. The chiral logarithms from each of the diagrams in fig. 1 are presented in table 2 
together with the total. 

Gasser and Leutwyler have calculated the chiral logarithms in the SU(3) theory as a 
function of q 2 [8]. In this case the power counting is different from that here, in that mx 
is also considered to be small. It is instructive to check our calculation by converting the 
SU(3) results to SU(2), using eq. (2.6) of ref. [8] and the expression for J in eq. (A. 7) of 
ref. [19]. Expanding the Gasser-Leutwyler results in powers of m 2 , we confirm that the 
total one-loop chiral logarithms in table 2 are indeed correct. 

From table 2 we now have all the ingredients to write down the NLO expression for 
the K — > tt matrix element at q 2 = 0. The expression is 



( 7r(P7r) I q^s | K(p K ) ) = F K p^ 



l-^L + c K ml 



1 - - L + On ml 



(3.7) 



where Fk = 1La\I j and F n = ILj^jj and c^ j7r are LECs. {Fk and F n should not be 
confused with the leptonic decay constants for which we use the notation fx and f n .) The 
chiral logarithm L is defined by 



2 / 2 

ml . / m 



and the dependence of L on /i is cancelled in expressions for physical quantities by the 
/i-dependence of the LECs (e.g. in eq. (3.7) the n dependence of L is cancelled by that of 
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Diagram 


Result 


%• 1(a) 
fig- 1(b) 
fig- 1(c) 
fig. 1(d) 


-^r 1 (PKfi -P-Kii) L 

2j f± {-hPK, - p^) + ^ l 


TOTAL 





Table 2: Tree level expression and the one- loop chiral logarithms for the K — > 7r matrix element 
at q 2 = 0. 



Cj,- and ex)- Eq. (3.7) implies that the chiral behaviour of the form factors is given by 

3 



/ u (0) = /+(0) = F + (1 - -L + c + mi) 
f-(0) = F^(l-^L + c.ml) 



(3.9) 
(3.10) 



where again F± and c± are LECs, given in terms of the parameters present in eq. (3.7) (for 
example, F± = \ (Fk ± F n ) ). 

Eq. (3.9) is the NLO SU(2) ChPT formula for extrapolating the lattice results for 
/°(0) = / + (0) which are obtained at unphysical values of the up and down quark masses 
to the physical point. The two LECs F + and c + need to be determined by fitting the mass 
dependence of the measured values of /°(0) to (3.9); the physical result of /°(0) is then 
readily obtained. Of course, using SU(3) ChPT the Ademollo-Gatto theorem [20] ensures 
that there are no LECs at one-loop order so that F + and c + are known and we can rewrite 
eq. (3.9) as 



/°(0) = /+(0) 



m 



K 



64vr 2 / 2 

,2 



1 + 



5- 12 log- 



4 Tfl T 

4 + 91og- + 31og^ 



(3.11) 



64vr 2 / 2 

The expressions for and c+ in Eq. (3.11) are valid only at linear order in m s ; the 
numerical results for / + (0) at small pion masses were found to lie below the one-loop 
SU(3) ChPT expression [3]. 

It is conventional for experimental results to be presented in terms of |K( S |/ + (0) and 
so we have concentrated above on the chiral behaviour of the form factors at q 2 = 0. We 
can perform a similar analysis for any value of q 2 with p n • px 
LECs, i.e. the F± and c± depend on q 2 . 



0(m 2 K ), but the effective 



f 

1 max 



4. Ke3 form factor at q 2 = gj 

We now turn our attention to the form factor / u (<Z max ), where q 



2 

max 



(rriK — fn n ) 2 . The 
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TOTAL 


^WO* (1 - ^L) + ^fp^ (1 - |L) + ^p^2L 



Table 3: Tree level expression and the one-loop chiral logarithms for the K — ► n matrix element 

& t 5max- 



tree- level diagram for the K — > 7r decay is drawn in fig. 1(a) and its contribution to the 
amplitude is given in the first row of table 3. By setting fx = 4 for example and recalling 
that 2L A1 = /W (see eq. (2.18)) we see that in the chiral limit f°(<7m ax ) = f (K) /f an d 
hence establish (1.4). 

In the remainder of this section we try to understand why the lattice results for / ((?max) 
in table 1 are significantly different from the value in the SU(2) chiral limit, f^/f, and 
seem to be approaching this value very slowly, if at all. At q^ iax the momentum of the 
external pion is small (px ■ Pit = raxm^) and so the counting of contributions in terms 
of powers of is simpler than at q 2 = 0. However, close to the SU(2) chiral limit 
(m n = 0) the corrections to the K — > it matrix element of the vector current are linear 
in m n (and not quadratic). To see this, consider for example, the matrix element in 
eq. (3.1), (ir(p 7T ) \ (D^D^K)^ D v (£ — £/)h \ K{j>k) ), which is now manifestly linear in m^. 
The coefficient of the linear term is not calculable directly in SU(2) ChPT. 

Below we study the chiral behaviour of /°(9max) i n three stages as follows: 

i) We start in section 4.1 by calculating the one-loop chiral logarithms, i.e. the correc- 
tions of 0(m\ log^) 2 ). This can be done within SU(2) ChPT. 

ii) In order to estimate the remaining terms we use SU(3) ChPT. In the second stage, in 
section 4.2 below we estimate the the coefficient of the linear term in by converting 
the one- loop SU(3) expressions to SU(2). 

iii) Finally, in section 4.3 we also estimate the quadratic terms using SU(3) ChPT . 

We will find that at g^ax the chiral corrections are very large and provide only a qualitative 
or perhaps a semi-quantitative, explanation of the observed chiral behaviour. Neverthe- 
less, the calculations confirm that the differences of the lattice results from f^/f are 
reasonable. 
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Figure 2: Sketch ofl — ll/4Lasa function of the mass of the pion. The three vertical lines 
correspond (from left to right) to the physical pion mass and to the lightest two masses in the 
simulation of ref . [3], 329 and 416 MeV respectively. \i was chosen to be m p = 0.77 GeV. 



4.1 The chiral logarithms 

The chiral logarithms from each of the diagrams in fig. 1 are presented in table 3. From the 
table, choosing the Lorentz index [x = 4 and neglecting terms of 0{m^ x ), with or without 
logarithms, we deduce that the chiral behaviour of the form factor is of the form: 



f0 ( 2 \ = ff_ 
J VSfmax/ j 



4 4vr/ (4vr/) 2 



(4.1) 



where Ai^ are low energy constants which depend on the strange quark mass but not on the 
light quark masses. Again one can readily verify that the coefficient of the chiral logarithm 
in (4.1) is indeed the result obtained by converting the general SU(3) formulae of Gasser 
and Leutwyler [8,19] to SU(2). 

The coefficient —11/4 is large (for example, at q 2 = in eq. (3.7) the coefficient of L 
is —3/4) and the term with the chiral logarithm does give a sizeable contribution in the 
region of pion masses between the physical one and that where the lattice simulations of 
ref. [3] were performed. In fig. 2 we sketch 1 — 11/4 L with the physical mass of the />meson 
as the scale fi and with / = 115 MeV which is the central value found in [9] . The sign of 
the chiral logarithm however, is such as to make the form factor decrease as the mass of 
the pion is decreased towards the chiral limit, which is the opposite of what is required to 
account for the difference between the measured values of /°(o , max) m table 1 and f^/f. 
Thus the chiral logarithms approximately double the size of the effect which should be 
explained. 

4.2 Linear term in m n 

We cannot evaluate Ai using SU(2) ChPT alone. To estimate whether the linear term in 
m-n in (4.1) can account for the difference of the measured form factors from f^/f we 
convert the SU(3) results of Gasser and Leutwyler [8, 19] to SU(2) ChPT. In this way we 
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Figure 3: Sketch of the expression in parentheses in eq. (4.2) with A2 = as a function of the mass 
of the pion (solid curve). The three vertical lines correspond (from left to right) to the physical 
pion mass and to the lightest two masses in the simulation of ref. [3], 329 and 416 MeV respectively. 
H was chosen to be m p = 0.77 GeV and / = 115 MeV. The dashed line represents the expression in 
eq. (4.3). 



can obtain an approximate value of Ax, using which we rewrite (4.1) as: 



2 ) 

max ' 



f 
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rnxm^ ( 14 20 . 4 8 r 
(4wU + Y l0g 3-9^ arCtan ^ 



V / 
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2m„ 
m K 



+ 



^2 2 , 

77V + 



(4tt/) 



(4.2) 



where rhx is the mass of the kaon in the SU(2) chiral limit. Eq. (4.2) represents an 
approximation for Ai since, within SU(2) ChPT, Ai contains higher powers of m s , whereas 
in (4.2) we have kept only those from one-loop SU(3) ChPT. Setting A2 = and neglecting 
higher order terms, we plot the expression in square parentheses in (4.2) as a function of 
the pion mass as the solid curve in fig. 3, where we have set f {K) /f = 1.28, / = 115 MeV 



and /i 



rrir 



We notice that the linear term in does indeed change the sign, the value 



of the form- factor does increase as we approach the chiral limit. The effect is too large 
however, and since the 0{m- K ) term is as large as 50-80% in the region where we have data, 
the stability of the chiral expansion is likely to be questionable. 

We write the converted expression from SU(3) ChPT as eq. (4.2) because this is the 
natural form for SU(2) ChPT. To illustrate that the result may depend significantly on the 
higher order terms we also plot, as the dashed curve in fig. 3, the expression 

'/<*> 



1 + 



f(K) 



11 
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14 20 , 4 

1 log - 

3 9 6 3 



9 



arctan 



f 



1 



2m,r 
fh K 



(4.3) 



which is equivalent to that in square parentheses (with A2 = 0) in eq. (4.2) at one-loop 
order in SU(3) ChPT but differs by terms which are powers of ffix/i^f)- We make this 
choice because it is the form obtained directly from one-loop SU(3) ChPT. The difference 
in the curves in the region where we have lattice data is about 25-30%, confirming that the 
uncertainties due to higher order terms are indeed likely to be large. 
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Figure 4: The curve is a sketch of /°(<7max) from eq. (4.2) with A2 = as a function of the mass 
of the pion. The three vertical lines correspond (from left to right) to the physical pion mass and 
to the lightest two masses in the simulation of ref. [3], 329 and 416 MeV respectively. For the curve 
\i was chosen to be m p = 0.77 GeV and / = 115 MeV. The black points are the lattice values from 
table 1 and the red points were obtained using SU(3) ChPT as described in the text. 



4.3 SU(3) ChPT 

Finally we use the Gasser-Leutwyler SU(3) ChPT results to estimate the effect of the 
chiral extrapolation of /^°^ from the quark masses used in the simulation in ref. [9]. In 
this case, at one-loop order, there is one LEC, the Gasser-Leutwyler coefficient L5, which 
is also the LEC which governs the SU(3) chiral behaviour of the ratio fx/fir] we can 
therefore use knowledge of the mass dependence of fx/fir (e.g. from lattice simulations) to 
evaluate the term proportional to L\. There is then no dependence on the scale \i. We use 
eq. (2.6) of ref. [8] to estimate the one-loop effects, but since the coefficients are large, the 
results depend on the precise procedure employed and on the choice of parameters (e.g. the 
physical value of f n or that in the SU(2) or SU(3) chiral limits), even though the differences 
are formally of higher order. As in section 4.2, we find that one-loop SU(3) ChPT predicts 
that the mass dependence of /°(gmax) * s steeper than that expected from table 1; the 
dependence on the pion mass is however, much less steep than in section 4.2. For example, 
using the physical value of the decay constant f n as the expansion parameter in the chiral 
expansion and taking the ratio of physical decay constants to be fx/fir — 1-20, we find 
that the one- loop prediction for the form factor at the physical quark masses is about 1.06 
(in the limit m n = the result, of course, is f^ K '/f — 1.28). Using the measured values 
of masses and decay constants from ref. [9], we find that, based on one-loop SU(3) ChPT, 
we would expect /°(gmax) to be about 0.94 for = 329 MeV and 0.90 for = 416 MeV 
as compared to 1.02 and 1.01 in table 1. To illustrate the flattening we plot in fig. 4, the 
expected chiral behaviour of f°(q^a ax ) using the chiral logs and linear term (as in fig. 3 but 
now multiplied by f^/f)- This is the curve in fig. 4. We also exhibit the predicted values 
from SU(3) ChPT as calculated above (the 3 red points) and the two lattice points from 
table 1. 

Thus one-loop SU(3) ChPT, with the procedure we have employed, can provide a semi- 
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quantitative explanation of the chiral behaviour observed in lattice simulations. Given the 
large effects we are finding at one-loop order, this is satisfying. We do stress however, that 
because of the large one-loop effects, the predictions are not very stable against varying 
the inputs into the chiral predictions, e.g. whether one uses the computed values of the 
decay constants (as we did above) or the values in the chiral limit. We have seen however, 
that the differences of the values in table 1 from f^/f are not unreasonable. 



5. Comparison with semileptonic B and D decays 

The chiral behaviour of semileptonic B — > ir decay amplitudes near g max was studied 
in refs. [11—14]. Chiral loop corrections were evaluated in standard chiral perturbation 
theory in [15] and extended to quenched and partially-quenched cases in [16,17]. For B- 
decays there is an additional scale, mt,, the mass of the 6-quark which is taken to be much 
larger than the typical hadronic scale Aqcd- Indeed the calculations in refs. [15-17] are 
performed by first taking the limit — ► oo (i.e. treating the 6-quark as being static) 
and then considering the chiral behaviour. In that case additional diagrams to those in 
fig. 1 have to be evaluated; in particular diagrams with B* propagators lead to one-loop 
contributions with chiral logarithms. This is because in the static limit the B and B* 
mesons are degenerate and so an on-shell .B-meson can emit a soft pion and the resulting 
S*-meson is also close to its mass-shell. The K and K* mesons on the other hand, are not 
degenerate and so the corresponding contributions are already contained in the diagrams 
of fig. 1 and the LECs. Thus the approach to the chiral limit depends on the order in 
which one takes the limits m& — > oo and m n — > and we return to this point below. 

The symmetry arguments used in section 4 which equate the semileptonic form factor 
/°( ( /max) m tne chiral limit to f^/f can be generalised to other flavours and in particular 
to semileptonic D and S-decays, so that, 



fD-tniOmax) 



and 



f° 

3 b- 



>7r(9max) 



/<*) 
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(5.1) 



where f( D ^ and are the D and l?-meson decay constants in the SU(2) chiral limit. 
This relation is also valid in the static limit (m& ~^ °°) [H — 13] and with l/m& corrections 
included in and f( B ) [14]. For a fixed finite value of m c or and for sufficiently 

small m w , the chiral corrections to the relations in (5.1) are given by eq. (4.1) with /W 
replaced by or f( B ) and the low energy constants Ai and A2 also depending on m c or 
m,b (as well as m s through strange sea-quark effects). In the static limit on the other hand, 
the approach to the chiral limit becomes [16] 
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(5.2) 



where gBB*w is the static BB*tt coupling. For fixed finite values of m&, eq. (4.1) may only 
represent the approach to the chiral limit at values of m n which are smaller than those 
accessible in lattice simulations and maybe even smaller than the physical value of m n . If 
that is the case then the behaviour of the form factor as a function of m n will have to be 
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studied either using the Heavy Quark Effective Theory (or NRQCD) as was done in the 
static limit in refs. [15, 16] or by paying explicit attention to the relative numerical values 
of Aq CD /m& and m w in QCD calculations. 

Eq. (5.1) provides an interesting check on the chiral extrapolations of lattice calcula- 
tions of leptonic decay constants and semileptonic form factors of B and D mesons. As 
these quantities are now being quoted with impressively small errors it is useful to have 
a constraint on the extrapolations. We postpone a detailed discussion of this issue to a 
future publication, but illustrate our point with an example. In fig. 14 of ref. [21] the 
authors plot the form factors fg^ n and as a function of q 2 in the chiral limit. Su- 

perimposed on the computed points is the fitted Ball-Zwicky parametrization [22] which 
suggests fB^n(lma,x) — 1-1> which is very considerably below the expected value of f^/f 
of greater than 1.7 or so 3 . This is an interesting puzzle which remains to be resolved. 
We have chosen this example because the f°(q 2 ) is helpfully presented in ref. [21] after 
the extrapolation to the chiral limit has been performed and so the results were relatively 
easy to interpret. A comparison of the chiral extrapolations in other studies with eq. (5.1) 
remains to be undertaken. 

6. Summary and conclusions 

In this paper we have studied the behaviour of the form factors as a function of the 
light quark masses (m u = ma) using SU(2) ChPT. At q 2 = 0, there is the subtlety that 
the final state pion is hard, nevertheless we have shown in section 3 that it is possible to 
calculate the SU(2) chiral logarithms. The one-loop expressions are given in eqs. (3.9) and 
(3.10). The coefficient of the chiral logarithm is small and we envisage that these formulae 
will be useful in extrapolating the results obtained for /°(0) in lattice simulations to the 
physical quark masses enabling a precise determination of the CKM matrix element V us . 

Following the procedure proposed in ref. [2] , lattice computations of K — > ir semilep- 
tonic form factors start with a very precise determination of / (<Zmax)- ^ n table 1 we see 
that the results at the values of the quark masses where the computations are performed 
are about 25% below the value in the chiral limit, f^/f- We investigated the chiral be- 
haviour of /°(9max) 111 an attempt to understand the difference of the lattice results from 
the value in the SU(2) chiral limit, f^/f. The coefficient of the chiral logarithms is of 
approximately the correct magnitude to account for this difference, but the sign is wrong; 
the 0(m 2 log(m^)) terms tend to make /°(<7max) decrease as we approach the chiral limit. 
We estimated the coefficient of the linear and quadratic terms in m, using SU(3) ChPT 
and found large effects, which lead to an increase in /°(9max) as m n decreases. In this way 
we obtain a semi-quantitative understanding of the difference of the lattice results from 
fW/f, but the large one-loop corrections prevent us from being able to determine the 
chiral behaviour with precision. 

There are a number of ways in which this exploratory investigation can be improved; 
one natural and necessary extension would be to perform the calculations at two-loop or- 

3 The points in the figure use the pre-erratum values from ref. [21]. However, the changes to the results 
for fs^-w are small enough not to affect the statement made here. 
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der in ChPT. The nature of the complementary relationship between the lattice and ChPT 
communities is changing in a very interesting way. Until recently, lattice computations 
were performed with quark masses which were at best marginally in a regime where ChPT 
could be applied (m^ > 500 MeV) and existing ChPT calculations were used to estimate 
the extrapolation to physical masses. Now as lattice calculations are being performed fur- 
ther into the chiral region (m n < 300 MeV), it is becoming possible to use the observed 
dependence on the momenta and particularly on the masses to determine the LECs and 
to test the range of validity and precision of ChPT (see ref. [9] for one such recent dis- 
cussion). We stress that in order for higher-order ChPT calculations to be useful in this, 
the expressions should be presented in terms of mass-independent LECs and with all the 
mass dependence exhibited explicitly. Of course, up to now, the primary aim of ChPT 
calculations has been to obtain predictions for physical quantities, i.e. for quantities at 
physical values of the quark masses, and for this it is sufficient to express the results with 
the numerical values for the physical decay constants and other quantities inserted into 
the expressions; for such a two- loop study of Kf 3 decays see ref. [23]. This prevents a 
determination of the full dependence on the pion masses and the calculated expressions 
cannot be used directly in conjunction with the lattice results. 

The calculations of the Kg^ form factors in ref. [3] were performed in unitary QCD, i.e. 
with valence and sea quark masses equal, and in this paper we studied the chiral behaviour 
using the standard unitary QCD. The use of partially quenched lattice simulations, in 
which the valence and sea quarks have different masses, is frequently a valuable tool in 
understanding the chiral behaviour of physical quantities. It would be a simple extension 
of the current work to evaluate the one-loop chiral logarithms in partially quenched QCD. 
Related to this is the use of partially twisted boundary conditions [24,25] in order to improve 
the momentum resolution in the calculations of form factors in general and enabling the 
evaluation of /°(0) directly without an extrapolation in q 2 [26]. The corresponding chiral 
and finite-volume corrections can also be evaluated. 

We wish to stress the generality of the relation between the semileptonic form factor 
fp^ n (Qma,x) an d the ratio of decay constants f^/f in the SU(2) chiral limit. It hold for 
all pseudoscalar mesons K, D and B and for sufficiently small values of the approach 
to the chiral limit at g^ax is given by the simple generalization of (4.1). Lattice calculations 
of the form-factors for semileptonic B —>■ ir decays are performed with the pion having a 
small momentum, i.e. at large values of q 2 . As usual, extrapolations in momenta and quark 
masses (as well as the lattice spacing) need to be performed and it is useful to use theoretical 
constraints to guide the extrapolations. One of these is the soft-pion relation in eq. (5.1) 
and its value in constraining the extrapolations remains to be investigated. Preliminary 
indications suggest that the conventional extrapolations lead to a value of /^-^^(^max) 
which is significantly smaller than f^/f, but this will be studied systematically elsewhere. 
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